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We show that dark matter annihilation around the time of recombination can lead to growing ion- 
ization fraction perturbations, that track the linear collapse of matter over-densities. This amplifies 
small scale cosmological perturbations to the free electron density by a significant amount compared 
to the usual acoustic oscillations. Electron density perturbations distort the CMB, inducing sec- 
ondary non-gaussianity. We calculate the CMB bispectrum from recombination, that is marginally 
observable by Planck. Even though electron perturbations can be markedly boosted compared with 
the Standard Model prediction, the dark matter effect in the CMB bispectrum turns out to be small 
and will be difficult to disentangle from the standard model in the foreseeable future. 



I. INTRODUCTION 

Galaxy and galaxy cluster dynamics and cosmological observations seem to imply that five out 
of six parts in mass of all matter in the Universe is composed of dark matter (DM), that is not 
accounted for by the Standard Model of particles. The particle nature of DM is one of the most 
intriguing puzzles of our time. Many efforts are invested in trying to solve this puzzle at direct 
and indirect detection experiments. It is important to identify astrophysical and cosmological 
processes where the particle interactions of DM, rather than its gravitational pull alone, may 
be of relevance. In this paper, we discuss a cosmological observable where DM interactions 
can modify appreciably the Standard Model prediction: linear perturbations to the ionization 
fraction of hydrogen. 

To summarize our main findings: DM annihilation can significantly change the evolution of 
linear cosmological perturbations to the free electron density, at and after the last scattering 
epoch of the Cosmic Microwave Background (CMB). Consistent with all current constraints, the 
effect can be as large as an O(10) enhancement for perturbations on small scales. Of course, it 
is not enough for DM effects to be large. For us to learn about it, the effect must also be visible. 
In the current paper, we looked for observable imprints of the electron perturbations in CMB 
non-gaussiainity. To no avail: even though electron perturbations can be markedly boosted, 
the main boost occurs slightly after last scattering, and on scales below the Silk damping scale, 
and so the non-gaussianity signal is small. In the rest of this introduction we expand on our 
motivations and lay out the structure of the paper. 

As is well known, dark matter annihilation or decay could modify the ionization history of 
the universe, giving rise to excess Thomson scattering compared to the Standard Model predic- 
tion fit- Tl3l] - This extra scattering damps power in small scale CMB temperature anisotropies 
and adds power in polarization. For DM annihilation, CMB constraints apply to the pa- 
rameter combination {{av)/m x ), where {av) is the velocity- weighted annihilation cross sec- 
tion and m x is the DM mass. For simple thermal freezeout models, the cross section is fixed 
(ctu)^ 3 x 10 -26 cm 3 s _1 GeV _1 and current CMB constraints based on precision measurements 
[14l - (l6l ] become important for DM masses below and of order 10 GeV. Planck data expected in 
the near future will either provide a detection or tighten the constraints. 

All of the CMB constraints in the current literature apply to temperature and polarization two- 
point correlation functions, or power spectra. Existing analyses usually consider the homogeneous 
ionization history, where the free electron density is taken to be a function of time only, n e = 
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n e (t). In some analyses (e.g., recently, [l2l fl7|). late time (z < 30) annihilation in non-linear 
halos is included but only in terms of a (model-dependent) boost to the smooth component. 

In contrast, our interest in the current paper is with linear cosmological perturbations to 
the electron density, characterized by S e such that n e — > n e (t) (1 + 5 e (x, t)). Power spectra 
are insensitive to these fluctuations, because in two point functions they enter at fourth (4th) 
order in the primordial curvature perturbations £ ~ 10~ 5 . Thus they cannot compete with 
the leading quadratic contributions. The leading observable where S e may play a role is CMB 
non-gaussianity, in particular the three-point function or bispectrum. Many inflationary models 
predict a very small primordial bispectrum so the first non-zero contribution may be due to 
deviations from linear evolution, that can be described via second order cosmological perturbation 
theory A finite first order S e produces second-order CMB inhomogencitics, transforming into 
finite anisotropy bispectrum. 

The paper is organized as follows. In [JTI]we calculate electron density perturbations. First, in 
Qll Al we review the standard recombination model and show that DM can make a sizable impact 
on the homogeneous hydrogen ionization fraction x e after last scattering, that may in fact be 
dominated by DM annihilation. As is well known, an 0(1) relative correction here is allowed 
experimentally, because (i) the residual ionization x e after last scattering is small anyway, and 

(ii) in the temperature power spectrum, the best measured CMB observable, excess Thomson 
optical depth is partially degenerate with the normalization A s and tilt n s of the primordial 
curvature fluctuations. 

In mi Bl we move on to cosmological perturbations. We build on the analysis of [l8| and give a 
semi-analytical derivation of the perturbation S e that applies at high redshift z > 700, relevant for 
CMB studies. It was realized in [l9l [2pj that at early times electron density perturbations follow 
an amplified copy of the baryon acoustic oscillations, with amplification factor ~ 5 corresponding 
to ionization waves. Extending the analysis to include DM annihilation we find a growing, non- 
oscillating, ionization mode that tracks the DM perturbations. The main result of this paper 
is that on small scales, this growing mode can boost 8 e by more than an order of magnitude 
compared to the Standard Model prediction, with peak amplification right after last scattering. 
Our calculations generalize earlier work that focused on later times long after recombination, see 
e.g. M- 

In ElIIII wc consider non-gaussianity. Several analytical and numerical studies have shown that 
the bispectrum from recombination is relevant for Planck and should be accounted for when 
searching for primordial non-gaussianity [22T - [28| . The leading sources appear to be second order 
metric and first order electron perturbations, inducing second order radiation terms. Rcfs. [22l — 
|2~H found the bispectrum induced by S e may be marginally observable by Planck. An order of 
magnitude amplification by DM annihilation then looks naively quite promising; we therefore 
compute the bispectrum induced by 5 e . In doing so, we have found the current literature lacking, 
specifically when it comes to perturbations on small scales. Our treatment of this problem will 
be reported separately in [29j . 

Our analysis shows that unfortunately, DM annihilation has little impact on the recombination 
bispectrum. The main reasons for this are: (i) the amplification to electron perturbations peaks 
immediately after CMB last scattering - largely missing the visibility window and hitting the 
early Dark Ages, instead; (ii) the DM effect rises on small scales below the Silk damping scale; 

(iii) in general, short wave electron fluctuations cannot affect long wave photon modes. This 
reduces the effect on squeezed triangles, where much of the signal-to-noise for the bispectrum is. 

DM annihilation can affect the evolution of matter temperature perturbations during the 
cosmic dark ages. While we do not pursue this avenue here, a natural means to try and detect 
the effect in the future would be through observations of 21 cm absorption [2l|, l30l - [33 |. Our 
calculation of the temperature perturbations extends previous analyses by properly accounting 
for the early initial conditions from the time of recombination. 

We conclude in t hVl In App. [A] we discuss the implications of non-local energy deposition by 
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DM annihilation. 

Throughout this paper we work with the following fiducial WMAP 7-year [35[ cosmology: 
fl b h 2 = 0.0226, fl DM h 2 = 0.112, h = 0.704, fl K = 0, r = 0.087, A s = 2.16 x 10~ 9 , n s = 0.963, 
with k p = 0.05 Mpc" 1 . 



II. DARK MATTER ANNIHILATION EFFECTS IN THE RECOMBINATION 

HISTORY 

In this section we compute the cosmological electron density to first order in perturbation 
theory We show that DM annihilation can cause a growing ionization mode, beginning around 
the time of recombination. This growing mode can boost electron perturbations by an order of 
magnitude compared to the case without DM annihilation. We start by a brief review of the 
homogeneous calculation and then move on to the perturbation analysis. 



A. Homogeneous calculation 



We follow the standard Peebles three-level atom formalism [36H39J and neglect helium ioniza- 
tion. The homogeneous free electron density is found by solving an effective Boltzmann equation, 

—JL + 3Hn e = Q e , (1) 
at 

with the DM ionization rate I x included in the collision term 

Qe = (f3 H e- ei2/TR (n H - n e ) - a H n^ C H + I x . (2) 

Here, the type-B recombination and ionization coefficients, as (7m) and /3jy(Tf{), are given 
in [38l . r39l | : £12 = 10.2 eV denotes the first excitation energy of hydrogen; Tm is the kinetic 
matter temperature; and nu is the total number density of hydrogen, ionized and neutral. The 
factor Ch denotes the probability of an n = 2 hydrogen atom to relax to the ground state before 
being photoionized and without exciting an adjacent ground state atom. It is given by 

r _ 1 + K H A H (n H - n e ) 

Uh ~ l + K H (A H + p H )(n H -n e y ( ) 

where A# ss 8.3 Hz is the two-photon 2s —> Is transition rate and the L a redshifting rate, 
k h 1 {z) = (8nH(z)/\l), is described in HHH. 
The DM ionization term is given by 

Ix = ~~a C-ion ■ (4) 

The quantity (w/a 4 e^f) denotes the proper rate per unit volume at which energy from DM 
annihilation is absorbed in the plasma, measured in units of the hydrogen ionization energy 
€h = 13.6 eV; we will return to this quantity shortly. The factor Ci on encodes the fraction of 
the absorbed energy which goes to ionizing the plasma. We use a crude parametrization of the 
partitioning of the absorbed energy between direct ionization, atomic excitation, and heating [l| : 
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We stress that Eq. ((5]) is a rough estimate. A more careful account of the energy partitioning 
deserves further study [UHll, but is beyond the scope of this paper. 

We now discuss (u/a 4 ) , the proper rate per unit volume at which energy from DM annihilation 
is absorbed in the plasma. For later convenience, where possible we will switch to work with 
conformal time drj = (dt/a). Thus here and in what follows an over-dot represents derivative 
with respect to r\. We will also use comoving coordinates ax ~ ^proper- Note however that unless 
stated otherwise, we still use proper particle densities, e.g. n e denotes the free electron density 
per unit proper volume, etc. 

Jumping ahead of ourselves for a minute by including spatial inhomogeneity, the comoving 
power density injected into the plasma is given by 1 



nj{x,r)) = a 4 (r))^-p 2 {x ,■!]). (6) 

The energy absorption rate is, in general, different from the injected power. At the epoch of 
interest, namely during and after recombination, particles coming out of an annihilation event 
can propagate over non-negligible time before their energy is absorbed by the plasma. The 
propagation time and distance depend on the particle type, initial energy, and time of injection, 
with final states of relevance being photons, electrons, and protons (with neutrinos trivially 
escaping indefinitely). The local rate of energy absorption is then given by folding the injection 
rate with some distribution, J 7 , specifying the propagation of the annihilation products, 

v r 

3„l 



{X,T))= I dlj d X f J r (x,3S , ,T),Tj')u in j(sS , ,Tj'). (7) 



We will have more to say about the quantity J- in App. For now, returning to the homo- 
geneous calculation, the spatial integral in Eq. © goes away and one is left with 

rn 

11(77)=/ dri'ui nj (r)')f dep (r),r)'), (8) 



where 



f dep (v,v') = —^ = I dV-F(af, a*, »?,»/) (9) 



-mj 



dlj 



describes the amount of energy absorbed by the plasma at time 77 per interval drj, after injecting 
initial energy €i„j at the annihilation time 77'. The time integral in Eq. ((5]) can be factored out 

as HiEllia, 

u(v) =Uinj(v)f(v)' (10) 
For time-independent velocity-weighted annihilation cross section, 

In App. [X] we discuss f(rf) in the context of concrete particle physics model examples. These 
examples illustrate the sensitivity of f{rj) to model details. Nevertheless, it is a reasonable 



For concreteness, we assumed Majorana DM and neglected possible time dependence in (uv). We remind 
that p x here still refers to the proper - not comoving - DM density. In the homogeneous limit it is given by 
a 3 Px(v) = a' 3 p x (ri'). 



5 




FIG. 1: Left: Ionization fraction, without DM annihilation (blue, solid; lowest curve) and with 
((av}/m x ) — 1.5 x 10~ 27 cm 3 /s/GeV (green, solid; middle curve) and {(av)/m x ) = 1.5 x 
1CP 26 cm 3 /s/GeV (red, solid; uppermost curve). We omit late time reionization in the plot. Dashed 
lines denote the floor approximation, Eq. ()13|) . Right: Relevant time scales; recombination rate per 
electron a.un e (blue solid thin curve), standard photoionization rate per electron f3n e _E12//Tfi (x~ L — 1) 
(black solid thick curve), DM ionization rate per electron I x / n e (red dot-dashed), and expansion 3H 
(green dashed). We take ((av)/m x ) = 1.5 x 10~ 26 cm 3 /s/GeV. 



approximation to take f(rj) as constant, /, over the time scale of recombination; for garden 
variety standard model final states, / ranges between 0.3 to 1 @, H3- Thus in this paper, as a 
rule, we simply absorb / into the definition of (av). 

Sample numerical solutions of Eq. (JTJ are depicted by solid lines in the left panel of Fig. [TJ 
where we plot the ionization fraction, 

Xe = — , (12) 

n H 

vs. rcdshift. As is well known, DM of mass m x = O(GeV) and annihilation cross section 
compatible with thermal freeze-out can have a significant effect on the ionization fraction after 
recombination. 

To understand the relevant processes, it is instructive to inspect the time scales appearing in 
Eq. (HJ. In the right panel of Fig. [1] we consider DM with ((av)/m x ) — 1.5 x 10 -26 cm 3 /s/GeV. 
We represent the DM ionization rate per electron by (I x /n e ), the standard ionization term (due 
to CMB photons) by /3jye _ei2 / Ti? '{x~ x — 1), and recombination by anrie- Notably, beginning at 
z ~ 1000 and down to z ~ 200, both DM ionization and recombination arc faster than Hubble 
expansion. Thus, for large enough annihilation rate, n e follows a quasi-equilibrium solution 
balancing recombination off DM ionization alone, 



a ^ = fe tej^M, (13) 

p b y 27 m x e H a H 

where we expanded to zcroth order in x e . This is the "floor solution" pointed out in (42j. 
Because an oc z~ 2 ^ 3 at the relevant redshift, a constant velocity weighted annihilation cross 
section results in 



1-3 



i/3 / (av) \ 1/2 / m x \-i/2 



1000/ V 3 x 10-26 cm3 / s J VGeV 



(14) 
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In the left panel of Fig. [JJwc depict the floor solutions by dashed lines. 

We caution the reader that DM with mass m x ~ 1 GeV and thermal relic annihilation 
cross section, as depicted by the red curves in Fig. [TJ is already excluded by CMB constraints. 
The strongest constraints we are currently aware of combine WMAP7 and SPT data to de- 
rive the bound {{av)/m x ) < 1.6 X 1CT 27 cnA^GeV- 1 at 95%CL WMAP7+ACT gives 
((av)/m x ) < 2.1 x lfT 27 cm 3 s- 1 GeV" 1 at 95%CL Q. WMAP7 alone gives ((av)/m x ) < 
5.7 x 10~ 27 cm 3 s" 1 GeV~ 1 at 95%CL [Io[. We thus use the example in Fig. |TJ merely to highlight 
the physics. 

The important point to take home is that DM annihilation can easily dominate the frac- 
tional ionization immediately after recombination. This remains true as long as the DM-induced 
ionization rate is comparable to the expansion right after recombination, namely, as long as 
((<jv)/m x ) > 10~ 27 cm 3 s _1 GeV _1 , corresponding roughly to the green curve in Fig. [TJ Some- 
what surprisingly, the current CMB constraints do allow for sufficient annihilation power. This 
is because most of the effect on the CMB temperature power spectrum on small angular scales, 
I > 100, is contained by an overall suppression factor C; — > e~ 2Ar Ci, where At denotes the 
excess optical depth due to the extra ionization. This overall factor is degenerate with adjusting 
the amplitude of the primordial curvature power spectrum, A s — > e 2Ar A s [42j . The amplitude 
degeneracy is not complete, and is ameliorated by polarization data; nevertheless, additional 
degeneracy with the primordial tilt n s and with other cosmological parameters leads to the fact 
that CMB constraints still allow a much larger role for DM annihilation at recombination than 
could naively be guessed. 

The main simple result of this paper can be understood directly from Eq. (fT3|). Generalizing 
to include cosmological perturbations, Eq. (|13|) tells us that in the quasi-equilibrium limit, we 
may expect the electron density perturbations to track DM perturbations, 



where the DM density contrast is given by 5 X = (Sp x /p x ). During and soon after recombination, 
DM perturbations on small scales are orders of magnitude larger than the corresponding baryon 
perturbations, 5 X 3> £&, because the latter are trapped in acoustic oscillations until the end of 
the baryon drag epoch while the former simply collapse gravitationally since horizon entry. Thus 
we may expect a large enhancement in the free electron density S e compared to the Standard 
Model prediction. Here we neglected several factors, including e.g. photon and kinetic matter 
temperature perturbations. However, in the next section we show that the simple reasoning 
behind Eq. (fT5"|) . motivated by the floor solution Eq. (fT5)) . is essentially correct. 

Finally, DM annihilation affects also the kinetic matter temperature, though around the time 
of recombination this effect is much less important than the effect on the ionized fraction. The 
equation for matter temperature is given by 2 




5n e 



S 



(15) 



dT M 



+ 2HT M 



7 



7 = 7c +7x> 



(16) 



dt 



with 



8aTdR 



n e T R (T R -T M ), 



(17) 



7c 



3m e c 




(18) 



2 We omit negligible corrections associated with photorecombination/ionization cooling/heating l43l . We thank 
Yacine Ali-Hai'moud for a discussion on this point. 
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where cir is the radiation constant. 

At redshifts z > 200, Compton scattering dominates and the matter temperature Tm tracks 
the CMB temperature Tr. At those early times, DM annihilation has negligible effect on the 
matter temperature. This is different than what we have just seen for the ionized fraction, and 
will carry over to the perturbation analysis in the next section. The reason is that, compared to 
the power available from DM annihilation, the CMB energy reservoir is intense but cool. DM 
annihilation can dominate ionization, because ionization can only feed off the deep Boltzmann tail 
of the CMB spectrum. The matter kinetic energy, however, is driven by Thomson scattering off 
the bulk of the CMB spectrum and thus DM looses this battle by a large margin. Quantitatively, 
we readily see this by inspecting the ratio of heating rates, 

1C _ 12 cr T c m x T R p HPl ^ f T M \ 
7 X (av) m H m e p 2 , e \ T R ) 

This comparison means that around recombination, when z ~ 10 3 and x e ~ 10 _1 — 10 -2 , DM 
annihilation cannot break the relation Tm = Tr, enforced by Thomson scattering. Only later 
at z < 10 2 and with x e ~ 10~ 3 — 10~ 4 , can DM annihilation compete with CMB heating. At 
this later time baryons kinctically decouple from the CMB and DM annihilation can change the 
evolution of Tm appreciably. 

B. Cosmological perturbations: inhomogeneous recombination 

We now compute the first order perturbations to the free electron density. Our aim is to refine 
the rough analysis leading to Eq. (fT5|) . As we discuss in more detail in Sec. HID what motivates us 
in pursuing this calculation, is that electron density perturbations during recombination induce 
apparent non-gaussianity in the CMB anisotropies as measured by a late time observer [l8|, [22| — 
|24| . This non-gaussianity can in principle be measurable by Planck and future experiments, 
offering a potential means to test our scenario. 

For simplicity we assume that the energy from DM annihilation is instantaneously absorbed 
by the plasma. As mentioned earlier, this can be a poor approximation, the extent to which it 
applies depending on model details. The smearing of energy absorption by the plasma leads to 
damping of small scale power. We analyze this issue in some detail in App. IA1 

We work in synchronous gauge, 

ds 2 = -a 2 (drj 2 - + h %3 ) dx l dx j ) , (20) 

fixing the gauge as usual by eliminating the DM velocity perturbations. Considering only scalar 
perturbations, we denote the trace and the trace-less parts of the scalar mode of hij by h and k. 
(These correspond to h and rj in the notation of Ma and Bertchinger [H|.) Our normalization 
for the primordial curvature perturbation is such that ^ = 1 on superhorizon scales. 

The Boltzmann and Einstein equations for metric (h, n), radiation (St r ), dark matter (S x ) and 
baryon density and velocity perturbations (8b and Vb), are not coupled to the electron and matter 
kinetic temperature perturbations (8 e and St m ) at first order. Therefore, for all perturbations 
other than 8 e and 8t m , we may use the usual set of Boltzmann and Einstein equations, given 
e.g. in Ma and Bertchinger 

Given the solutions for h, K, 8t r , 8 x , 8b, and Vb - amounting to the usual transfer functions 
- we use them as sources for the linearized electron and matter temperature perturbations. 
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Starting with the results of [18|, we add DM annihilation to obtain: 



A 



• h 2ik Q 7 (sr- f dlog-f \ n e 6 e + n b 8 b \ 

°T M = - o ^~v b + ■ — — > — d x - Tm ; • (22) 

3 3 (n b + n e )T M \^ \ologX J n b + n e j 



Here X = {H,n e ,n b ,n X7 Tju,TR}, and 3 5h = — [S b /3aHJ denotes the perturbation to the 

baryon velocity divergence, as measured by a comoving local observer. 

It is straightforward to solve Eqs. (j21j|22l) numerically. However, a simplification occurs if one 
suffices with computing observable effects in the CMB 4 . As discussed in the previous section, 
the matter temperature is clipped to the radiation temperature around last scattering. This 
partially carries over to the perturbations: 5t m ~ St r , all the way until the end of the baryon 
drag epoch when St m rises by compression as the baryons fall into the DM potential wells. By 
the time 5t m finally breaks loose of 5t r , then, the Thomson optical depth for photons is small 
and the electron perturbation has little residual effect on the observed CMB anisotropy. 

Setting St m = St r , we can write a direct integral solution for S e , 

S e (k, rj) = J ' dr]'g e (k, rf) exp (- £ d V "T e (rf^j , (23) 

^ = aQe / _ fd\ogQe\\ g e = 6 b+ ^j:(^%)6 X ,, 

n e \ \d\ogn e J J n e ^ \d\ogX J 

where the sum over X' now does not include n e . The initial time ?yi n it is chosen early enough 
so that 5 e (k, r]i n it) = S b (k,r]i n it). Eq. (|23p allows us to obtain S e directly and quickly, reading all 
other perturbations from the numerical code CAMB [4^-[13l i it agrees well with the full numerical 
solution to Eq. ((2I~j) throughout and for a good while after recombination. 

As expected from the discussion in the previous section, the calculations confirm the presence 
of a growing ionization mode sitting on top of the usual baryon acoustic oscillations (B AOs) . We 
now examine this result in some detail. 

We begin with an eye towards observability in the CMB. In Fig. [2] we fix the wave number of 
the perturbation and examine its time dependence. Green line shows the electron perturbation 
with ((av)/m x ) = 3.75 x 1CP 27 cm 3 s _1 GeV _1 ; blue shows the electron perturbation with no DM 
annihilation. For reference, we show also the baryon and DM density perturbations in grey and 
black, respectively. Note that for 5 e and <5b, we plot the absolute value of the transfer functions, 
while the DM transfer function is positive. The grey shaded band depicts the full width half- 
maximum (FWHM) of the visibility function (taken here with no DM annihilation). In the left 
panel we fix k = 0.04 Mpc -1 , corresponding to observed anisotropy multipolc I ~ fcr/o ~ 600, 
where 770 ~ 1.4 x 10 4 Mpc is the conformal time today. In the right panel we fix k = 0.3 Mpc -1 , 
corresponding to I ~ 4200. 

The growing ionization mode due to DM annihilation is clearly visible. This mode grows 
towards, and finally catches up with the DM perturbations, eventually amplifying S e by more 



3 Note that the J? —dependence of Q e is contained in the Ch factor for La escape, Eq. J3). For the matter 
temperature, we have (d log ■y/d log H) = 0. 

4 The reasoning here may well need to be modified if one aims to address physics at later epochs, e.g. for 21 
cm analyses. For such late time effects one needs to solve Eqs. t'2]\i and (122 B simultaneously - as we do were 
required in this paper. 
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FIG. 2: Perturbations vs. conformal time. The curves are: S x (solid thick black); |<5f,| (dashed thin red); 
\S e \ with m x = 8 GeV (solid green); \S e \ with no annihilation (dashed thick blue). For the m x = 8 GeV 
curve, thermal relic annihilation cross section is assumed. Grey shaded band denotes full width half- 
maximum of the visibility function. Left: wave number k = 0.04 Mpc, corresponding to / ~ 600. Right: 
k = 0.3 Mpc, corresponding to I ~ 4.2 x 10 s . Synchronous gauge; £ = 1 on superhorizon scales. 



than an order of magnitude compared with the Standard Model prediction. However, it takes 
finite time for the quasi-cquilibrium configuration to manifest itself, particularly so on larger 
scales; this causes much of the amplification effect of 6 e to only take place after CMB last 
scattering. In addition, note that dragging the electron perturbation in the positive direction 
towards the DM perturbation, can actually lead to suppression of the magnitude - in absolute 
value - for perturbations that enter the last scattering surface with negative amplitude. This is 
seen in the left panel of Fig. [2] 

In Fig. [3] we study the scale dependence. In the left panel, we show the k dependence of the 
perturbations in a snapshot close to the peak of visibility r\ = 285 Mpc. On the right we focus 
on the half-maximum width, rj ~ 310 Mpc. The effect is larger for larger wavenumber, as DM 
perturbations on smaller scales enter the horizon earlier and have more time to collapse before 
recombination, leading to more efficient ionization. In addition, as noted above, the growing 
mode becomes significant only somewhat after the peak of the visibility. We caution the reader 
again that the red curve with m x = 2 GeV is excluded experimentally by CMB data, and is only 
shown here for illustration. 

In Fig. U] we show 5 e vs. redhsift at times significantly after CMB last scattering, again for 
two different wavenumbers. Fig. [4] concerns the deep Dark Ages; the detailed dynamics at last 
scattering, crucial for CMB analyses, is merely seen as small wriggles around z ~ 1100. Here 
we confirm the naive estimate of Eq. (|15[) . that says that for high enough annihilation power, 
the electron perturbation reaches quasi-equilibrium and roughly sticks to the DM perturbation. 
In this regime, S e is roughly independent of the DM mass and annihilation rate and is boosted 
by a factor 2-3 compared with the Standard Model prediction. Reducing the annihilation power 
below ((av)/m x ) ~ 10~ 27 cm 3 s~ 1 GeV _1 , as seen for m x = 36 GeV, causes ionization and 
recombination to drop below the Hubble rate, freezing the ionization fraction below the quasi- 
cquilibrium value. 

During the dark ages, the relevant future probe of DM annihilation would be in the absorption 
of 21 cm radiation [2ll. [3fj| - [33 |. (See Ref. [48| for a review.) Here, the relevant quantity is the 
matter temperature entering the computation of the spin temperature (49l . [5fj| . In the left panel 
of Fig. [5] we plot the matter temperature perturbation at z = 200 as a function of wavenumber. 
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peak visibility: r|=285 Mpc half-width visibility: T)=310 Mpc 
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FIG. 3: Electron perturbation vs. wave number, at peak (left) and half-maximum (right) of visi- 
bility. Thermal relic annihilation cross section is assumed. The curves are, from top to bottom: 
m x = 2, 8, 18 GeV, and no annihilation. Synchronous gauge; £ = 1 on superhorizon scales. We caution 
the reader that the red curve with m x = 2 GeV is excluded experimentally by CMB data, and is only 
shown here for illustration. 



k=0.05/Mpc k=0.3/Mpc 




FIG. 4: Electron perturbations vs. redshift. Left: wavenumber k — 0.05 Mpc, corresponding to / ~ 700. 
Right: k — 0.3 Mpc, corresponding to I ~ 4.2 x 10 3 . Thermal relic annihilation cross section is assumed. 
The curves are, from top to bottom: m x = 8, 18, 72 GeV, and no annihilation. Synchronous gauge; £ = 1 
on superhorizon scales. 

We learn that a factor of ~ 2 enhancement in St m can arise from DM annihilation. In the 
right panel, we plot the ratio of the baryon and matter temperature perturbations to the DM 
perturbation as a function of redshift, for fixed wavenumber k = 0.1 Mpc . As an aside, 
using Eqs. (|2HI22[) we can solve for the matter temperature perturbations from the correct initial 
conditions at recombination down to the deep dark ages. Doing this, we note that the scale- 
independent relation ST u (k,T]) = s(rf) Sb(k,rf) with s(r]) independent of k, assumed e.g. in [2l| 
and later references, is violated at 0(1). 

Previous analyses of CMB non-gaussianity induced by perturbations to the free electron density 
around recombination (l8l . [221424^ , have found a level of non-gaussianity that could be marginally 
detectable by Planck. We have seen that DM annihilation could boost small scale electron 
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FIG. 5: Left: matter temperature perturbation vs. wavenumber. Thermal relic annihilation cross 
section is assumed. The curves are, from top to bottom: m x = 8, 18, 72 GeV, and no annihilation. 
Right: baryon and temperature perturbations relative to dark matter perturbation. The top (dashed) 
curve shows 5b/8 x . The bottom four curves are as in the left panel. Synchronous gauge; £ = 1 on 
superhorizon scales. 



perturbations by a sizable amount. It is thus interesting to assess the bispectrum when DM 
annihilation is taken into account. 

As we show in Sec. IIII1 despite the amplified electron perturbations, the DM annihilation 
effect on the CMB bispectrum is small. This comes about from three unfortunate reasons. First, 
photon diffusion acts to erase power on small angular scales, where the DM effect is pronounced. 
Second, as we saw recombination has a finite response time. By the time the DM-induced 
amplification reaches its full swing, photon last scattering is mostly over. And third, Thomson 
scattering cannot transmit power from a short wave electron perturbation down to a long wave 
photon anisotropy 5 . Thus, short wave electron perturbations do not contribute directly to the 
bispectrum in the squeezed limit, where much of the signal-to-noisc is contained. At the end of 
the day, our results indicate that the large boost to 5 e will be very difficult to detect in the CMB 
even if the recombination bispectrum is measured, at least in the foreseeable future. 



III. NON GAUSSIANITY: CAN WE OBSERVE ENHANCED SMALL SCALE 
ELECTRON PERTURBATIONS IN THE CMB? 

Assuming gaussian initial conditions, any observed non-gaussianity and, in particular, a finite 
bispectrum comes about at second order in perturbation theory. In this section we write ap- 
proximate formulae for the second order temperature anisotropics in the presence of (first order) 
electron density perturbations, and estimate the induced bispectrum. 

We stress that the analysis presented here is meant as a rough estimate of the observabil- 
ity of the effect highlighted in the previous section. We defer a more comprehensive analysis 
to [29]. There, special care is given to electron perturbations on small scales, that have not been 
accounted for by existing analytical studies. The results then apply gencrically and no special 



5 In the relevant limit of photon number conservation. 
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treatment is required to include DM annihilation, beyond utilizing the modified 5 e computed 
above. 

To save the reader from disappointment: our final answer to the question posed in the title 
of this section is negative. We find that even an O(10) enhancement, compared to the Standard 
Model, in small scale electron perturbations, will leave only a very subtle imprint on the bispec- 
trum. As most of the DM annihilation effect during last scattering is concentrated on such small 
scales, the CMB bispectrum will likely not provide means of detection. 



A. The bispectrum 

Our notation for the homogeneous (unperturbed) differential and integrated optical depth and 
the visibility function are given by 

rno 

f(n) = -aca T n e , r(r/) = dr)'t(r)'), 5(77) = -e^ T(v) T(r]). (24) 



We write the Fourier space temperature anisotropy as 

Q(k,r],h) = 6 (1) (fc,?7,n) + e (2) (fc,?y,n). (25) 

We neglect second order metric perturbations. Then, the first and second order anisotropics 
today are given by the line of sight (LOS) solutions: 



6 (1) (fc,?7o,n) = / d77e ifc ^ ( ^ o) S (1) (fc,7?,n), (26) 
Jo 

e^(k, Vo ,n) = j d V e ik ^-^ g( V ) (S Sg (k, V ,n) + S^(k, V ,h)j , (27) 
where we define fik = k ■ n. The source terms are: 

S^(k,r],n) =.g(^e[ ) 1) (fc) + /i fc 4 1) (fc)-i7' 2 0i fc )n (1 )(fc) + 2a(fc)^ + e~ T (k(k) + S(*)) + ga(k), (28) 

s Sg (k, v ,n) = J g) (eS 1 '® + ft^ 1 '® - ^^(^inf^?) - ewip)) , (29) 

(k, n,n) =e 2 \k) + h- vf ] (k)-\v 2 {iik)rt 2 Xk). (30) 

We define a = (h + Qk)/2k 2 (5lJ and suppress the r\ dependence on the RHS for clarity. The first 

order baryon velocity perturbation is assumed to be irrotational, v^\k, rj) = kvjp (k,rf). 

In Eq. ([50]) we neglect vector and tensor contributions (m = ±1,±2, respectively; note that 
the rotational velocity vanishes in the bispectrum). Then, in the source terms, we can use the 
same Legendre multipolc decomposition for first and second order terms, 

9j(M) = ^ / dnVi(n-k)Q(k,r),ti), (31) 

where Vi(x) are Legendre polynomials and with similar decomposition for the polarization Op;, 
feeding into LI = 2 + Opo + ©P2- Note that for the second order perturbation, ©( 2 \ the 
Legendre decomposition of Eq. (|31|) contains only part of the information because at second 
order, azimuthal symmetry around the wave vector no longer holds. Thus when we compute 
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spherical harmonic coefficients, we will need to use the full Yi m transform for <d^ 2 ' (see Eq. (f3"3")l 
below). Nevertheless, focusing on scalar contributions, the Lcgcndre moments in Eq. ([3"U)) suffice 
to compute the second order source because there azimuthal averaging occurs by Thomson 
scattering. 

Physically, the contribution S$ g comes about from perturbing f and t in the first order so- 
lution, Eq. (|26p . This term is equivalent to perturbing the visibility function, up to corrections 
proportional to the integrated Sachs- Wolfe effect that are only relevant on large scales. As written 
in Eq. (|2l?)) . this term is ready to deploy in an explicit calculation of the three-point function. 

The contribution contains the second order feedback. Namely, it includes the actual effect 
of the electron perturbation on the photon field, rather than, as in the previous term, the effect of 
perturbing the way we see that field today. We defer a derivation of to subsequent work [2§] . 
Our results for this term extend and improve the analysis of (23), and disagree with [24|. 

In what follows, we derive the bispectrum contribution due to the visibility term S$ g and 
discuss the effect of DM annihilation. We then discuss qualitatively the contribution due to the 
second term S^ 2 K 

The spherical harmonic coefficients a/ m = + aP^ are given by 



42 =4nJ ^H)'e, (1) (£;, 



d 3 k 

^ / dhYC m {h)& 2 \k,h lVo ). (33) 



(2*) 

Using Eqs. (|32II33[) we compute the bispectrum, Bf^^ m = (ae imi ae 2m2 ae 3Tn3 ) . After the dust 
settles, we find that the contribution due to Sg g leads to the following result 6 : 

B^ltlm 3 = Gt[t 2 m 3 X — / ^(l) (fii MgiM + five permutations) , (34) 

n Jo 

gt( V ) = J dkk 2 P{k)Q ( p (k,r) ) ji[Km - v)Mk,v), 

Mv) = (-i)' / dkk 2 P(k)e£\ k , m ) £(2i' + i)(2i" + i) J +v+l " 3 v[Km-v)] 



( 9^(k,r,)-e { ;\k, v ) nW(M) n r U1 x , n (i )( , A 

x I <Vi — hd ; »2 (l-dpo)(l-<5rx)0iV(K,»7) I • 

In Eq. (|34p . by 8 e (k, 77), 0;(fc, 77), etc. we refer to transfer functions, namely, we mean to have mod 
out the random initial curvature perturbation £z from S e (k, r/), Qi(k, rj). The variables 9^ = ikv^ 
and 6 1 = 3fcOi are as in (44|. We assume gaussian adiabatic initial curvature perturbations, 
with power spectrum 

(^) = (2Tr)H^(k + p)p(k). (35) 

Finally, the gaunt coefficient is 



, hh i 3 = / (2^ + l)(2£ 2 + l)(24 + l) (iiiih) h h h 
' mim2 " 13 « 4^ [ o o I 1 mx m 2 m 3 



Our derivation of Eq. (1341) follows closely that of [23I. |24H . but the result disagrees with theirs in a few terms. 
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FIG. 6: The function gi(r)) of Eq. (JMJ), capturing the DM annihilation effect in the visibility contribution 
to the recombination bispectrum. We plot gi (in absolute value, units of fiK) vs. I, with a snapshot 
at peak visibility (left) and half-maximum visibility (right), for different DM masses assuming thermal 
freeze out cross section. 



The information about the electron density perturbation is contained in the function gi(rf). 
DM annihilation affects gi(rj) at high I. In Fig. [6] we plot gi(rj) (in absolute value, fiK units) 
vs. /, with a snapshot at peak visibility (left) and half-maximum visibility (right). We examine 
different DM masses, assuming standard thermal freezeout annihilation cross section. As it turns 
out, the growing ionization mode develops too late to show up significantly during last scattering 
where the visibility is large, and on angular scales that are too small to be accessible to current 
and upcoming experiments. Only in the 2 GeV case, that is already excluded experimentally, a 
noticeable effect occurs at peak visibility and on scales I < 2000. Imposing existing constraints, 
the effect is small even on scales I > 2000, beyond the reach of Planck. 

Obviously the effect of DM annihilation in the bispectrum is small, on scales I < 2000 where 
other secondaries such as point sources are under control (see e.g. [52[). Nevertheless, for com- 
pleteness we estimate the signal-to-noise. The signal-to-noise ratio for the detection of the 
bispectrum, for an experiment with detector noise Nt and covering a fraction of the sky / s ky, is 
given by 



(37) 



with A^^., = 1,2,6 for zero, two and three equal £'s. The reduced bispectrum is obtained by 
summing over m-modes, 

5 W3 = E h (38) 

' — ' \ mi mi ) 

The S/N arising from Eq. is plotted in Fig. [71 for the case without annihilation and for 
thermal freezeout annihilation cross section with m x = 8 GeV. Where S/N is not negligibly 
small, the annihilation scenario is essentially indistinguishable from the Standard Model. In 



computing S/N, we follow [22j and only include angular scales I > l m in — 100, that are well 
within the horizon during recombination. This procedure is meant to eliminate the sensitivity of 
the result to second order metric perturbations, that have not been included in our estimate of 
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the bispectrum. We checked that using smaller values of l m in has negligible effect on the results. 
We consider cosmic variance limited (CVL) experiment with zero detector noise and / s k y = 1 . 



10° 
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FIG. 7: Signal-to-noise ratio for a CVL experiment. Blue - no DM annihilation, green - m x = 8 GeV 
with thermal freezeout annihilation cross section. The plot includes the bispectrum from perturbed 
visibility only, Eq. ((34)) . 



Lastly, so far we analyzed the bispectrum contribution coming from perturbed visibility, S$ g , 
but ignored the contribution of the second order source S^ 2 \ This is not because the contribution 
due to S^ 2 * 1 is small; it is in fact comparable to the Ss g term that we discussed [22j]. To our 
knowledge, a correct analytical estimate for the bispectrum contribution from is yet to be 
published. This analysis is motivated regardless of the specific implications for DM and we take 
it up in [29j|, where we also explain where current estimates [22I. |24| are lacking. However, even 
without the detailed answer for the contribution due to S^ 2 ' , there is a simple physical argument 
which makes clear that DM annihilation can only slightly modify this contribution from the 
Standard Model result. To understand this argument, note that S^ 2 ' encodes the cumulative 
effect of the electron perturbation on the photon multipoles up until last scattering. 

Now, a rather accurate description of CMB anisotropics can be obtained within the tight 
coupling approximation [53[ , where the effect of Thomson scattering on the CMB multipoles is 
packaged into effective diffusion (or Silk) damping, 

e 0l i(M)~6 0l i(M)e ^- ( 39 ) 

Here, Oo,i(fc,f?) are the photon monopole and dipole, obtained from the Boltzmann equation 
deleting Thomson scattering and neglecting all other multipoles, and 



k 2 D ( V ) J d71 2r{l5 ' 1 + R 



R 2 



(40) 



is the diffusion scale with R = (4pb/3p 7 ) and c~ 2 = 3(1 + R). 

Ref. [22j used Eq. (|3T))) to derive a rough estimate of the bispectrum contribution due to 
S^ 2 \ limiting the analysis to electron perturbations 5 e (k) on very large scales, k <C kjj. By 



7 To compare with expected performance of Planck, there the beam size will cut the growth of S/N above 
I ~ 1500. 
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construction, their derivation can not strictly apply to small scale S e and thus can not be used to 
assess quantitatively the effect of DM annihilation. Nevertheless, Ref. [22| pointed out that small 
scale S e has negligible contribution to the bispectrum, and this observation remains qualitatively 
correct. The reason to this is simple: electron perturbations on scales smaller than the diffusion 
mean free path, can not affect diffusion damping. Thus Eq. (|39[) . by packing Thomson scattering 
into an effective diffusion coefficient, already implies that short wave electron perturbations 
cannot lead to big effects. 

We comment that additional contributions associated with S^ 2 ' exist, that are not captured 
by diffusion damping; these contributions are identified in 29] in terms of perturbations to the 
photon-baryon sound speed and baryon drag. However, these terms are less significant than the 
diffusion effect and do not change the results appreciably. 

IV. CONCLUSIONS 

We compute linear perturbations to the free electron density 6 e , including the effect of dark 
matter (DM) annihilation. We find a growing, non-oscillating, ionization mode that tracks the 
DM perturbations. The main result of this paper is that on small scales, this growing mode can 
boost 8 e by more than an order of magnitude compared to the Standard Model prediction, with 
peak amplification right after last scattering. The kinetic matter temperature is also affected 
with potentially 0(1) corrections from the Standard Model prediction during the cosmic dark 
ages, relevant for 21 cm observations. 

CMB power spectra are insensitive to these linear electron density fluctuations. The lead- 
ing observable where S e may play a role is CMB non-gaussianity, in particular the three-point 
function or bispectrum. There, a first order electron perturbation feeds into second order, non- 
gaussian temperature multipolcs. Several analytical and numerical studies have shown that the 
bispectrum from recombination is relevant for Planck and should be accounted for when search- 
ing for primordial non-gaussianity. Refs. [2214241 ] found the bispectrum induced by S e may be 
marginally observable by Planck. An order of magnitude amplification by DM annihilation then 
looks naively quite promising; we thus computed the bispectrum induced by S e . In doing so, we 
have found the current literature lacking, specifically when it comes to perturbations on small 
scales. Our treatment of this problem will be reported separately in [29j . 

We find that the non-gaussianity signal is small, very difficult to disentangle from the Standard 
Model by any current or upcoming experiment. This is because even though electron perturba- 
tions can be markedly boosted, the main boost occurs slightly after last scattering and on scales 
below the Silk damping scale. 

While the prospects for observation in the CMB look slim, our analysis does show that signif- 
icant 0(1) changes to the ionization history of the Universe may be caused by DM interactions 
during the early cosmic dark ages. In particular, an 0(1) enhancement of electron density and 
matter temperature perturbations, with power rising on small scales similarly to DM pertur- 
bations, would follow from DM annihilation. Similar conclusions were found for later epochs 
relating to DM halos; our analysis extends these findings to the early linear regime. A natural 
observational tool to try and detect these effects in the future is 21 cm radiation. 
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Appendix A: Non-local energy deposition 

Not all of the energy injected by DM annihilation is absorbed by the plasma, and the deposition 
of the part that is absorbed takes non negligible time. The eventual energy deposition occurs 
when the annihilation shower becomes an electromagnetic cascade, with electrons and photons 
cooling down to the ~ keV range where ionization and heating takes over. Here we describe this 
effect in the homogeneous limit and then proceed to estimate the implications when cosmological 
perturbations are included. We also make contact with computations of [lTj to illustrate the effect 
in some concrete model examples. 

Ref. [TtJ computed the object T(z,z'), defined separately for different initial energy e, n j and 
done for electrons and photons: 

T(z,z')dz = — §i. (Al) 

This relates to our conventions via 

f dep (r,,r,')=T(z,z')H(z). (A2) 

In the homogeneous limit, the effect of non-local energy absorption is then encoded by the 
function f(rf) of Eq. {TTTJ), 

f(r,) = jT dtf (a/a'f f dep (77, r,') = J dz |M f I±iQ * T (z, z >). (A3) 

The case of DM decay, or of time-dependent (av) , is a simple generalization of Eqs. (TTT]) and (|A3|) . 

For the purpose of computing cosmological perturbations, both the time and the spatial smear- 
ing of the energy deposition are relevant. Linearizing Eq. ([7J and moving to Fourier space for 
clarity, 

5u d ep{k,rj) = Uinjiv) J dr)'(a/a') 2 {sf dep {k,r],r]'^ +2 5 x (k,r)')J 7 (k,r),r)') S ) . (A4) 

For the homogeneous part of we used F (x, x + f, 77, 77') = .F(0, 7^,77,77'), with Fourier trans- 
form T (fc, 77, 77'). For the perturbation in J 7 , it is natural to generalize the quantity fdep, 

J d 3 x'5T (x, x , 77, 77') = 6 f dep (x, 77, 77') , (A5) 

with Fourier transform <5/rf ep (fc, 77, 7/'). 

The term 5 fdep comes from various non-DM perturbations. For instance, electron density 
perturbations affect the cooling time of energetic photons in the electromagnetic shower following 
DM annihilation. At late times, z < few hundreds, baryonic perturbations are as large as DM 
density perturbations. Then, we expect the 5 fdep term to be as relevant as the DM 6 X term 8 . In 



This can be relevant e.g. for 21 cm analyses. 
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this paper, however, we restrict our interest to the recombination epoch where, for modes inside 
the horizon, the DM perturbations S x are much larger than all other baryonic (and metric) 
perturbations. To obtain basic understanding of the physics, it is safe to neglect the Sfdep 
term in Eq. (|A4[) . In addition, again around recombination and for modes inside the horizon 

a'S x (k,r/) « aS x (k,r)'). Using these observations we can write, 

rn 

6u dep (k,r]) w 2u inj (T])5 x (k,T)) drj '{a/a')? {k,r],r)') . (A6) 

Jo 

To proceed further, we need information about the model dependent distribution T . Let us 
consider simple examples. 

• Instanteneous deposition: consider DM with mass m x ~ 100 MeV annihilating to e + e _ . 
Close to the time of recombination at z ~ 10 3 , the electrons cool quickly by inverse Comp- 
ton (IC) scattering on CMB photons, with a comoving cooling scale 

k c » 2.5 • 10 3 (z/W 3 ) 3 (e/GeV) Mpc" 1 . (A7) 

For ~ 100 MeV electrons, this gives k c ~ 260 Mpc -1 corresponding today to angular 
resolution I ~ 10 6 , beyond our current ambition. Thus these electrons quickly and locally 
deliver their energy to IC photons with typical energy e 7 ~ Je^CMB ~ keV. The ionizing 
photons quickly deposit their energy in the plasma. For this model it is a reasonable 
approximation to assume instantaneous deposition, 

T (£, x 1 , 77, rf) « 5^ (x - of )6( V - rf), (A8) 

leading to 5u dep (x,r]) « 2u in j(rj)8 x {x,rj). 

In the left panel of Fig. [S] we explore fdepiViV') f° r the 100 MeV, XX e + e~ example 
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Deposition smearing: consider now DM annihilating to e + e~, but with larger DM mass 
m x ~ GeV. The initial electrons still cool quickly by IC scattering, however, some fraction 
of the energy will now go to IC photons with energy e 7 ~ MeV. These MeV photons 
are non-ionizing; they must cascade down by Thomson scattering to the keV range before 
they can be absorbed by the plasma. Thus, some fraction fi ns t of the initial annihilation 
energy will be deposited locally, but the remaining 1 — fi ns t will be smeared over significant 
distance and time. Consider an ansatz for the deposition smearing, 

.F(f,f,77,V) w hnsM)^ ) (x-x l )5{r 1 -,i) (A9) 

3 

1 



+ (1 - finstW)) 1 e 2 ^"""' ) /rfe P , 7 (W) 

Here, fdep,-y describes the energy loss rate of the secondary IC photons in the plasma. As 
the photons of interest have MeV energy - comparable, but not much exceeding the self 
energy of electrons in the plasma - the spatial smearing here should be quite similar to 
CMB diffusion damping. In particular, around the time of recombination we can estimate 



3^70 k 2 D ■ 



9 We thank Tracy Slatyer for providing us with high resolution grids of her results of energy deposition. 
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FIG. 8: Left: energy deposition for XX ~~ ^ e e > m x = 100 MeV, at z = 1080. Right: same, but for 
m x = 1 GeV. 



For the non-local part of the energy deposition, then, we have diffusion damping. Plugging 
into Eq. 

/ rv _2k 2 \ 

6u d ep(k,r)) ^ 2u inj (r/)5 x (k,r)) \f inst (ri) + J dr)'(a/a') (1 - /i ns t(»/)) /dep, 7 (v, v') e *°J(All) 

In the right panel of Fig. [5] we plot the function fdepiViV') f° r the 1 GeV, XX e + e~ 
example. A certain fraction of the injection energy can be attributed to a narrow peak 
immediately attached to the annihilation time. However, much of the absorbed energy 
exhibits significant deposition time, much larger than for the previous example and relevant 
in comparison with the time scale of recombination. This extended deposition can be shown 
to arise from the shallow behavior of fdep,y 

To summarize, annihilation energy deposition on small scales is damped by photon diffusion. 
Considering Eq. (|A11[) . neglecting the time dependence of all factors compared with that of 
/dep. 7 {VtV') m the T]' integral, we can estimate 

Sude P {k,T)) ^2u inj {t])S x (k,r]) (finstiv) + (1 - finstiv)) Mv) e~ 2fc2/fc °) , (A12) 

where / 7 (?y) = J drj' (a/a/) fdep,-/ {ViV')- I n analyzing DM annihilation as a source for cosmologi- 
cal ionization and matter temperature perturbations, we should thus keep in mind that a model 
dependent, but potentially no n- negligible fraction of the annihilation power in DM density per- 
turbations on small scales, k > kjj ~ 0.15 Mpc , is washed out and does not source ionization 
or temperature perturbations on these scales. 

Finally it is clear that, analyzing distributions such as in Fig. [5J the quantities fi ns t and 
fdep,-y (or more generally, the smeared component) can be readily extracted. Once this is done, 
Eq. (jAllj) or (|A12[) can be used to calculate the energy absorption damping effect for cosmological 
perturbations. In this paper, due to the model dependence of the processes involved, we will not 
go into these details. 
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